Mechanochemical action of the dynamin protein 
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Dynamin is a ubiquitous GTPase that tubulates lipid bilayers and is implicated in many membrane 
severing processes in eukaryotic cells. Setting the grounds for a better understanding of this biolog- 
ical function, we develop a generalized hydrodynamics description of the conformational change of 
large dynamin-membrane tubes taking into account GTP consumption as a free energy source. On 
observable time scales, dissipation is dominated by an effective dynamin/membrane friction and the 
deformation field of the tube has a simple diffusive behavior, which could be tested experimentally. 
A more involved, semi-microscopic model yields complete predictions for the dynamics of the tube 
and possibly accounts for contradictory experimental results concerning its change of conformation 
as well as for plectonemic supercoiling. 
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I. INTRODUCTION 



In order to operate efficiently, living cells must con- 
stantly maintain concentration gradients of various chem- 
ical species and isolate some of their components. One of 
the many different biological processes required to main- 
tain this traffic is membrane fission, by which a cell mem- 
brane compartment is split into two or more topologically 
distinct parts. A fundamental protein involved in most 
membrane fission events is dynamin, which has been pro- 
posed to be a "universal membrane fission protein" [1]. 
Dynamin and its analogues are found in cellular processes 
as diverse as clathrin-coated endocytosis, phagocytosis, 
mitochondria and chloroplasts reorganization, cell divi- 
sion and virus immunization in organisms ranging from 
mammals to yeast and plants [I] 15] . In most of these pro- 
cesses, two separating membrane compartments end up 
being linked by a thin membrane neck which is difficult 
to sever, since the membrane must be strongly curved 
before it can be pinched off. Dynamin-like proteins lo- 
calize at this neck and participate in its breaking, thus 
completing the fission. Although mutations in such an 
important protein are often lethal, defective dynamin 
or dynamin analogues have been shown to be involved 
in human diseases such as the optical atrophy type 1, 
the Charcot-Marie- Tooth disease and the dominant cen- 
tronuclear myopathy O HI |5] . 

The role of dynamin in tube fission was first suggested 
by results showing the importance of its Drosophila 
analogue in endocytosis ^|. Dynamin is recruited by 
clathrin-coated vesicles, possibly through membrane- 
mediated elastic interactions [7], and self-assembles into 
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short (a few helical repeats) helical constructs on the 
cell membrane necks localized at their base [8 . Much 
longer (thousands of helical repeats) helical dynamin 
polymers have been observed in cell-free environments, 
either wrapped around microtubule templates [9j, or in 
solution |10[ lllj . Purified dynamin also polymerizes 
around negatively-charged lipid bilayers [lll[T3], deform- 
ing liposomes into dynamin-coated nanotubes, simply 
termed "tubes" in the following. Electron micrographs 
suggest that these tubes are hollow, i.e. filled with water 
[HI [T5j [16] . Dynamin is a GTPase and therefore cataly- 
ses the hydrolysis of guanosine triphosphate (GTP) into 
guanosine diphosphate (GDP) and inorganic phosphate 
(Pi). This highly exoenergetic reaction (~ 25fcBT per 
GTP molecule in a typical cellular environment) is simi- 
lar to the hydrolysis of adenosine triphosphate (ATP), 
which fuels most known molecular motors and many 
other cellular processes [17], and has been shown to be 
necessary for endocytosis [IS]. Self-assembly of dynamin 
has been linked to a dramatic increase of its GTPase ac- 
tivity [mUnilll], thus suggesting that GTP hydrolysis 
by self-assembled dynamin drives membrane tube fission 
during endocytosis [TJ |5J [55] . 

During the past decade, experimental evidence indi- 
cating that tube breaking involves a mechanochemical 
action of dynamin has been accumulated. Electron mi- 
croscopy has shown that the geometry of the dynamin he- 
lical coat changes upon GTP hydrolysis. However, there 
is still a controversy as to whether both the pitch and ra- 
dius of the helix shrink [HI [H [TSJ (111 [SS] (see Fig. [I]; a)) 
or the pitch increases at constant radius [24l [25] . In the 
absence of any other protein than dynamin, this change 
of conformation is sufficient to drive tube breakage when 
the end points of the tube are attached to a substrate 
[13] . However, no fission of freely floating tubes is ob- 
served [53]. More recently, optical microscopy has been 



used to investigate the dynamics of the tube's conforma- 
tional change and breaking [26j . In these experiments, 
dynamin-coated membrane nanotubes are grown from a 
lamellar phase of a suitable lipid mixture. Then GTP is 
injected in the experimental chamber (typically in a few 
tenths of seconds) . The initially rather floppy tubes then 
straighten, revealing a build up of their tension. If one 
end of the tube is free to fluctuate, this tension results in 
the retraction of the tube. If both ends are attached, the 
tube breaks. If polystyrene beads (diameter 260-320 nm) 
are attached to the dynamin coat, rotation is observed af- 
ter GTP injection, showing that GTP hydrolysis induces 
not only tension but also torques in the tubes. The typ- 
ical time scales involved in the rotation of the beads and 
the breaking of the tubes are roughly 3 s after GTP in- 
jection. 

In this article, we present a theoretical model for the 
relaxation of long dynamin-coated membrane nanotubes 
accounting for the above-mentioned experimental results. 
We beheve that a quantitative description of the tube dy- 
namics will help to understand the mechanism by which 
dynamin severs membrane tubes. This is a much-debated 
question for which several models have been proposed 
[22j . Since little quantitative information about the mi- 
croscopic details of the dynamin helix is available, we 
choose a coarse-grained (hydrodynamic) approach. In 
this framework, we do not need to speculate about the 
unknown microscopic details of the non-equilibrium be- 
havior of the tube: its dynamics is characterized by a few 
phcnomcnological transport coefficients. 

In the next three sections, we present the building 
blocks of our formalism by decreasing order of general- 
ity. In Sec. [n] wc consider only the symmetries of the 
system and write the most general hydrodynamic the- 
ory compatible with these symmetries. We then argue 



in Sec. Ill that one hydrodynamic mode is much slower 
than the others. This leads to simplifled equations de- 
scribing this mode. In Sec. |IV[ we present two micro- 
scopic models of the equilibrium properties of the tube 
aimed at describing two possible experimental situations. 
This allows us to solve the equations of motion and make 
predictions about the tube dynamics. In Sec.lV] we com- 
pare these predictions to experimental results and thus 
justify some of our assumptions. A tentative account of 
the differences in the conformational changes of dynamin 
reported in Refs. [131 EH HSl HI US] on the one hand and 
Refs. [231 [2S] on the other hand is also given. Finally, we 
discuss the generality of our model and its implications 
for membrane nanotube flssion in Sec. I VII 



II. HYDRODYNAMIC THEORY 

In this section, we derive equations of motion for dy- 
namin/membrane tubes based on the symmetries of the 
system. We also restrict our study to the long length 
and time scales, thus constructing a hydrodynamic the- 
ory. More speciflcally, we focus on the so-called hydro- 
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Figure 1: Schematics representing the geometry of the tube. 

(a) The tube comprises two fluids h and m, here pictured in 
difi'erent colors. It is invariant under a rotation around z by 
an arbitrary angle 6 followed by a translation by p6 along 
z (we refer to this property as helical symmetry with pitch 
2ttp in the text). Moreover, the system is assumed to be 
invariant under a rotation of vr around the x-axis (the sys- 
tem is non-polar) . The latter transformation is equivalent to 
a reversal of polar coordinates {9,z) — > {—6,—z). The nu- 
merical values measured in Ref. [15 for the radii and pitches 
in the relaxed/constricted states are indicated on the figure. 

(b) Representation of the same helix in the 9, z plane of cylin- 
drical coordinates with periodic boundary conditions on 6. 
The system clearly has a broken translational symmetry in 
the direction of the translucent arrow. It is described by a 
broken-symmetry variable u^e obeying Eq. ( |10| , which can 
be understood as a conservation law for the number of stripes 
visible here. The associated reactive current is Vh/p — i^h, the 
projection of the velocity of the stripes on the direction of the 
arrow. 



dynamic modes, which are spatially inhomogeneous ex- 
citations of the system away from equilibrium with the 
following properties [27] : 

1 . the amplitude of these excitations are small enough 
for the system to remain weakly out of equilibrium 
in the sense of Ref. [28 , 

2. the wave vector q and the pulsation uj{q) charac- 
terizing the spatial inhomogeneity of the hydrody- 
namic mode are such that 



lim a;(g) = 0. 

9—0 



(1) 



The q -^ limit corresponds to excitations over length 
scales much larger than the microscopic length scales of 
the system. Typically, we consider inverse wave vectors 
of the order of the tube's length: q~^ ~ several 10 /im. 
This is indeed much larger than the typical microscopic 
length: the tube radius r ~ 10 nm. The hydrodynamic 
theory thus involves a coarse-graining of the system at 
the scale of a few tens of nanometers, and so the tube 
must be treated as a one-dimensional object. 



Let us now state the hypotheses underlying our hy- 
drodynaniic theory. We consider a one-dimensional sys- 
tem comprising two fluids which we refer to as fluids h 
(representing the helix) and m (the lipid membrane). In 
agreement with electron microscopy data [T31[T3[Tn], we 
assume that the system has a helical symmetry with pitch 
27rp and is non-polar, as shown on Fig. ^fa). 

In the following, we identify the relevant variables de- 
scribing this system and derive an expression for its en- 
tropy production. Introducing an active term represent- 
ing the input of free energy in the form of GTP, we write 
the constitutive (flux/force) equations for the tube. To- 
gether with conservation laws these equations eventually 
yield the hydrodynamic modes of the system. 



A. Conservation laws and hydrodynamic variables 

The first step in building a hydrodynamic theory relies 
on conservation laws. 

We assume that no exchange of membrane or dynamin 
occur with the aqueous medium surrounding the tube on 
the time scale of the change of conformation ^26J . There- 
fore, the masses of fluids h and m obey the following 
conservation equations: 



with proportionality coefficients {7ij}j 



dtPh 
dtPm 



-'^{PhVh) 
-^{PmVm), 



(2a) 
(2b) 



where V is the differentiation operator with respect to z. 
Vfi and Vm are the velocities of fluids h and m respectively 
and ph and pm their mass densities (masses per unit of 
z length). We now define the mass fraction of ft, as $ = 
Ph/p, the mass density of the whole tube p — ph + Pm, 
the linear momentum density of the tube g = phVu + 
PmVm = pv with V the center-of-mass velocity and the 
diffusion flux of h relative to the center of mass of the 
tube J = ph{vh — v)- The conservation laws expressed in 
Eqs. (l2| can be re-written as 



dtp 
dt<i> 



-V.g 



-^VJ, 



(3a) 
(3b) 



It is shown further below that the inverse relaxation times 
of p and <i> go to zero with vanishing q, meaning that p 
and <& are hydrodynamic variables. 

Let I be the angular momentum density of the tube. 
The conservation laws for g and I are the force and torque 
balance equations. There are two contributions to the 
force (resp. torque) applied to a tube element: first, the 
divergence of a (resp. r), the linear (resp. angular) in- 
ternal stress of the tube; and second, the external force 
(resp. torque) due to the coupling of the helix dynamics 
with the hydrodynamic flow that it induces in the sur- 
rounding aqueous medium. For simplicity, we model this 
"friction against water" as a force and a torque linearly 
dependent of v^ and the angular velocity flh of fluid h 



dtg = 


-VfT- 


- IzzVh - 


- IzB^h 


(4a) 


dtl = 


-Vr- 


- lezVh - 


-jee^h- 


(4b) 



Since Vh = {g + J/^)/p, we can replace Vh hy g + J/^ 
in Eqs. Q after rescaling the friction coefficients in the 
following way: jzz -^ Plzz and 7^2 -^ Plez- Since the 
friction between the two fluids is a local phenomenon 
(which does not depend on the wave vector g), the quan- 
tity f2/i — f2„j relaxes to zero in a non-hydrodynamic time 
(that does not respect Eq. (fTl)). When studying hydro- 
dynamic time scales, we can thus replace fi^ by ri„i or 
equivalently by 51 = 1/ 1 with / the tube's density of mo- 
ment of inertia. Redefining 7^9 — > I^zq and 70^ -^ Ijee, 
we replace flh hy I in Eqs. Q. 

It is now apparent in Eqs. Q that in the presence of 
friction against water, g and / relax to the solutions of 
the following equations: 



7z^ ( 5 + ^j^ j + Izei = -V(T 
7^2 ( 3 + ^ j + leel = -Vr 



(5a) 
(5b) 



over times of order l/jtj- Since the jijS do not depend on 
q, these times are not hydrodynamic times (they do not 
go to infinity when q vanishes). The linear and angular 
momentum densities g and I are therefore not hydrody- 
namic variables and are given by Eqs. (|5|, which are the 
force and torque balance equations in the overdamped 
regime. 

The conservation of energy reads: 

dte + Vf = V (-jzz ig + ^j " Izel 

+n (-^ez (g + ^j - leelj , (6) 

where e is the energy density and j^ the current of energy. 
The right-hand side accounts for energy dissipation by 
friction. 

Although no other conservation law than Eqs. (l3|, ((5| 
and ([6| exist in the system, its hydrodynamic descrip- 
tion is still incomplete. Indeed, the system has a broken 
continuous symmetry similar to that of smectic-A liquid 
crystal phases (Fig. [ijb)). Just as in the case of hquid 
crystals [23], we use a strain tensor component Uze to 
describe this symmetry breaking. We define it in the 
following way: let 9{z,t) be the angular displacement of 
the intersection of the helix with the plane located at al- 
titude z at time t (therefore is a eulerian coordinate). 
The strain is then defined as: 



Uze{z,t) = — (z,i). 



(7) 



As any broken-symmetry variable, Uz9 obeys a relation 
similar to a conservation law. To show this, we first note 



that in a fully reversible situation 

dteiz,t)^nh~-. (8) 

p 

Differentiating this equation with respect to z, one finds 

'vh 



P 



(9) 



where Vh/p — ^h is the reactive current of Uz9 (see also 
Fig. [iTb)). In the presence of dissipation, one must add 
a dissipative part X to this current [30 , hence 



dtu.e = -wx-w( — -nh 
\p 

^-WlX+^-±^-^). (10) 

pp I ; 



B. Entropy production 

The great simplicity of hydrodynamic theories can be 
tracked back to the fact that in the long-time limit, one 
locally describes the state of a potentially very complex 
system using only a few conserved quantities. Indeed, on 
time scales going to infinity with the size of the system, 
all the microscopic (fast) degrees of freedom have relaxed 
and the system is locally in a state of thermal equilibrium 
in the thermodynamic ensemble defined by the conserved 
quantities. 

Let us apply this idea to the tube in the general case 
where friction against water is not necessarily present (i.e. 
the 7yS can be zero, in which case g and I are hydrody- 
namic variables). The state of local thermal equilibrium 
is entirely characterized by the six quantities g, I, p, <&, 
Uzd and e. Equivalently, we can consider a homogeneous 
tube of length V and study it in the thermodynamic en- 
semble {P,L,M,^,Uz9,T,V), where P = Vg, L = VI 
and M = V p are the total linear momentum, angular 
momentum and mass of the tube; T is the temperature. 
The total differential of the free energy of the tube reads: 

dF = vdP+VL<:\L + pdM+M^iA<^ + Hduze-S<lT-pdV. 

(11) 
This equation defines the total and exchange chemical 

potentials p and /ig, the reactive stress h = H/V, the 

entropy S — Vs and the equilibrium pressure p. Note 

that the velocity v is thermodynamically conjugated to P 

and the angular velocity f2 to L. This implicitly assumes 

that the free energy is the sum of a kinetic energy and a 

static contribution: 



F = V 



pv 






+ F^{M,<^,Uze,T,V), (12) 



where the free energy of the tube in its rest frame can 
also be written as 

i^o(M,$,M,e,T,y) = F/o(p,$,u,e,T). (13) 



Table I: The fluxes, forces and signature of the forces un- 
der two symmetry operations: "time symmetry" denotes the 
time-reversal symmetry and "spatial symmetry" refers to the 
reversal of the polar coordinates (6*, z) -^ (— ^, —2) defined in 
Fig. [3a). 



Flux 


Force 


time 


symmetry 


spatial 


symmetry 


a-p-h/p 


Vv 




- 




+ 


T + h 


vn 




- 




-1- 


J 


V^e 




+ 




- 


X + J/(p$p) 


Vh 




+ 




- 


Q 


VT/T 




+ 




- 




Ap 




+ 




+ 



Using an extensivity argument, one shows that 

p = vg + ni + pp + Ts - e (14) 

and proves the local form of the fundamental equilibrium 
thermodynamic relation: 

Tds ~ ^vdg — ildl — pdp — ppe<i^ ~ hduzo + de. (15) 

Inserting the conservation equations Eqs. (I3|, Q, (l6| 
and ( 10 1 into (llS]) and using ( 14 1, one finds the following 



form for the local entropy production of the system 



T 



ds 
di 



= - 0--P- 



P 



Vv 



X 



{T+h)Vfl 
JVpe 

p<^pj 

vr 
- Q^, (16) 

where Q is the heat current in the z direction and where 
higher-order terms in the displacement from equilibrium 
have been dropped. 



C. Constitutive equations 



The right-hand side of Eq. ( 16 1 is the sum of five terms 



each of which is the product of a flux and a force as 
displayed in Table |lj These fiuxes and forces vanish at 
thermal equilibrium. Also, according to the second law 
of thermodynamics, entropy production is always posi- 
tive. Therefore, close to the equilibrium state, the fluxes 
depend linearly on the forces through a positive definite 
matrix. In addition to this positivity condition, other 
constraints exist on the relationships between fiuxes and 
forces: 

First, entropy production is invariant under spatial 
symmetry operations that leave the system unchanged. 
Therefore, fiuxes and forces of opposite signature un- 
der the spatial symmetry (0, z) — > {—9, —z) defined in 



Fig. [T|^a) cannot be coupled. This property is a spe- 
cial case of the Curie principle, which states that in an 
isotropic system, couplings between fluxes and forces of 
different tensorial characters are forbidden. In this con- 
text, the quantities displayed in Table |l] that are odd 
under the transformation {9,z) -^ {—d,—z) are analo- 
gous to vectors and those that are even are scalars or 
second-rank tensors [2H]. 

Time-reversal symmetry imposes another set of con- 
straints. Each flux can be written as the sum of a dis- 
sipative part, which has the same time symmetry as the 
conjugate force, and a reactive part with the opposite 
symmetry. Dissipative couplings occur between fluxes 
and forces having the same time-reversal symmetry. Con- 
versely, reactive couplings relate fluxes and forces of op- 
posite symmetries. According to Onsager's relations, the 
matrix of dissipative couplings is symmetric and the ma- 
trix of reactive couplings antisymmetric |31j . 

While deriving Eq. (16), we have made sure that its 



right-hand side involves only entropy production terms, 
and no entropy exchange or energetic effects. Therefore, 
the fluxes in this equation are dissipative and have a van- 
ishing reactive part. Taking into account the symmetry 
constraints discussed above, we obtain the following set 
of constitutive equations: 



CT — p — h/p = 

T + h = 

J = 

X + J/{p^p) = 



-■qz^v - aVn, (17a) 

-aVv - T]gVn, (17b) 

-AV^e - bVh ~ cVT/T, (17c) 

-bVne->^yh-dVT/T, (17d) 

-cV^e - dWh - kWT/T. (17e) 



The coefficients in front of the forces are so-called phe- 
nomenological transport coejjicients. They are a priori 
unknown coefficients that depend on the microscopic de- 
tails of the problem. 



D. Discussion of the phenomeno logical coefficients 

In the spirit of the present article, those phenomeno- 
logical coefficients that are relevant to the relaxation of 
the system should be determined experimentally. The 
only way to calculate them a priori would be to use a 
detailed microscopic model, which would require a bet- 
ter knowledge of dynamin than we have. However, in the 
next few paragraphs, we try to interpret the origin and 
give typical orders of magnitude of these phenomenolog- 
ical coefficients. 

The coefficient r/^ > can be identified as a length x 
surface viscosity, where "length" denotes a typical micro- 
scopic length of the tube, for instance its inner radius r ^ 
10 nm. Similarly, Tyg > is a (length)'^ x surface viscosity. 
Assuming that the effective characteristic surface vis- 
cosity of the tube is close to that of a lipid bilayer, 
namely of order 5 x 10~^kg.s~^ |31], we estimate that 
r]z ~ 10~^^kg.m.s~"^ and r]g ~ 10~^^kg.m^.s~^. 



The momentum transfer from translational to rota- 
tional degrees of freedom is described by a. This trans- 
fer is allowed since the tube is chiral. The amplitude 
of these effects is constrained by the positivity of the 
matrix of phenomenological coefficients, which imposes 

The coefficient A > relates a gradient of chemical 
potential to a diffusion flux. By analogy to Pick's law, 
we expect it to be proportional to a diffusion coefficient. 
To better interpret A, let us set all other phenomenolog- 
ical coefficients to zero. The hydrodynamic dissipation 
then comes only from the homogcnization of the helix 
mass fraction <i> at fixed mass density p, and therefore 
involves a relative flow between the two fluids h and m. 
In this scenario, the source of dissipation is obviously 
the friction between the two fluids. One can therefore 
interpret A as the inverse of a helix/membrane friction 
coefficient. In order to calculate an order of magnitude, 
we propose an oversimplified friction mechanism inspired 
by Ref. |33j . where it is shown that dynamin inserts 
into the outer leaflet of the membrane bilayer. In this 
naive model, we assume that the outer membrane mono- 
layer is attached to the helix and that the energy dissi- 
pation comes from the sliding of one monolayer against 
the other. Experimentally, one measures typical friction 
coefficients for the relative sliding of lipid monolayers of 
order f3 — 10* Pa/(m.s^^) J4J- Let us consider a motion- 
less isothermal (VT = 0) cylinder of membrane of length 
L surrounded by an undeformed (h — 0) helix of dynamin 
moving at velocity v^ under the influence of a chemical 
potential gradient difference pe = LV p^ between the ex- 
tremities of the cylinder. The mass flow of helix in such 
a system is phVu — P^Vh, hence the tube receives a net 
power V = pep^Vh fro m th e reservoirs located at each 
end of the cylinder. Eq. (17c) and J = p^{l — ^)vh entail 



V/L = p^^^{l — ^)vyX. Assuming that this power is 
entirely dissipated by the friction between the membrane 
and the helix implies V/L = 2i:r(3vf^ and eventually 



A = 



$2(1-$) p2 

VJ3 



27r 



1.1 X 10~2^kg.m~\ 



(18) 



where the values at equilibrium p^o = Po*&o — 3.7 x 
10^13 ]jgjn-i andp„o = Po(l-^o) ^ 3.8x lO'^kg.m"! 
are calculated from the molecular mass of dynamin [21] 
and the number of dynamin monomers per helix turn |15j 
on the one hand, and from the typical mass per unit area 
of a lipid bilayer [35, on the other hand. 

The coefficient A > has properties similar to those of 
A but only exists if the system has a broken-symmetry 
variable. If the system under study were a crystal, we 
would interpret this coefficient as related to the phe- 
nomenon of vacancy diffusion, i.e., the displacement of 
mass without change in the periodic lattice. In our sys- 
tem, unlike in a crystal, there can be two independent 
diffusion coefficients oc A and oc A even if the creation 
of vacancies (and therefore possibly the breaking of the 
helix) is forbidden. Indeed, one does not need to create 
holes in the helix to displace mass without disturbing the 



periodic lattice of Fig. [T|b) : this can be done by chang- 
ing the radius of the helix. As far as orders of magnitude 
are concerned, we only assume in the following that the 
transport phenomena associated with A are not much 
faster than the ones associated with A. 

In the following, we consider the system as isothermal. 
This condition can be enforced by making the thermal 
conductivity k > go to infinity, which implies that any 
thermal gradient relaxes instantaneously. In this k ~* oo 
limit, we can drop Eq. (17el as well as the last terms of 



Eqs. ([ITcf and (jlTdl. 

Eventually, b, c and d describe couplings between the 
three diffusion phenomena described above. Such cross- 
effects give rise for instance to the so-called Soret and 
Dufour effects. As in the case of a, the positivity of 
the matrix of phenomenological coefficients sets upper 
bounds on their values. 



E. Coupling of GTP hydrolysis or binding to the 
dynamics 

We have now developed a complete formalism for the 
dynamics of a passive, non-polar, diphasic helix submit- 
ted to external friction. However, the system considered 
in this article is not passive since nucleotide (i.e. GTP 
or GTP analogue in this context) hydrolysis by dynamin 
or at least binding to dynamin is required for conforma- 
tional change. In the following, we introduce this exter- 
nal free energy source using arguments similar to those of 
Ref. [36] . Instead of deriving a whole new formalism tak- 
ing into account the conservation of GTP, GDP and Pi 
and all the chemical reactions involving them, we model 
the presence of GTP in the experimental chamber by a 
spatially homogeneous "chemical force" A/x, where A/i 
stands for the free energy provided by the hydrolysis (or, 
arguably, binding) of one GTP molecule. 

From Table |I| we see that the spatial symmetry of A/i 
only allows it to couple to Eqs. (17a I and (17b). We 



also note that time-reversal symmetry imposes that these 
couplings are reactive. Neglecting thermal diffusion as 
discussed in Sec. |IID| we obtain a modified set of consti- 
tutive equations: 

fT - p - h/p = -77^Vw - aVn + ^^A/i, (19a) 
T+h = -aVt; - r/eVf} -I- ^eA/i, (19b) 



J = -AV/ie - b\7h, 



X 



p^p 



-bS/pe - AV/i. 



(19c) 
(19d) 



F. Hydrodynamic modes 

The hydrodynamic relaxation modes are studied by 
linearizing the equations of motion around the state of 
thermal equilibrium. By definition, all thermodynamic 



forces vanish at thermal equilibrium, and in particular 
A/i = 0. Let dp = p - po and i5<i> = <I> - <I>o be the de- 
viations of the mass density and of the mass fraction of 
fluid h from this state. Combining the conservation equa- 
tions of Sec. II A with the constitutive equations Eqs. ( 19 1 



yields dynamical equations relating Sp, u^e and 6^ with 
g, I and the reactive (equilibrium) forces p, h and p^. 
In the overdamped regime, g and / can be eliminated 
according to Eq. ([s]). Close to equilibrium, the reactive 
forces are linearly related to the state vector of the sys- 
tem X = {Sp,Uzg,S^) through a susceptibility matrix x- 



(20) 




Uz9 I = XX- 

(5$ 



It is therefore possible to write linearized dynamical 
equations for x in a closed form. In Fourier space and to 
leading order in the wave vector q, it reads 



iojx = -q^ (A'j-^B'' + aA XX 



(21) 



where A''', S'', 7 and A"^ are matrices. A^ and B"^ de- 
scribe reactive couplings, the elements of 7 are the ^ijS 
defined above and A"^ contains dissipative phenomenolog- 
ical coefficients. See Appendix [A] for details. According 
to Eq. (21 ), the system has three diffusive hydrodynamic 
modes. 



III. LONG TIMES DYNAMICS FOR THE 

HELIX/MEMBRANE FRICTION-LIMITED 

REGIME 

The results of Section |TT] allow a full description of 
the hydrodynamic behavior of the dynamin/mcmbrane 
tube. For instance, to predict the relaxation of a helix 
with some known initial and boundary conditions on the 
hydrodynamic variables {Sp^Uze^S^), one should diago- 
nalize the matrix 



M ■ 



U^'-f-^B' 



A' 



X 



(22) 



and solve three diffusion equations along the directions 
defined by its eigenvectors. We denote the eigenvalues 
of M as Di > D2 > -D3. Unfortunately, this diago- 
nalization yields very lengthy expressions from which no 
intuitive picture of the dynamics can be deduced. Nev- 
ertheless, we show in this section that all experimentally 
observable features of the tube's dynamics can be faith- 
fully described by more convenient simplified equations. 
The matrix M is the sum of two terms. This re- 
flects the fact that the dynamics of the tube involves 
two sources of damping: A^'y^^B'^x describes the fric- 
tion against the outer water and A'^x is associated with 
dissipation mechanisms internal to the tube, such as he- 
lix/membrane friction. We now make an estimate of the 
orders of magnitude of these two effects. Assuming that 



the friction of the tube against water is that of a rigid 
rod of radius re (defined as the external radius of the 
dynamin coat) and length L yields [571 



Iz 



27rr/p 



In 



-0.72 



Ize = lez 



0, 7ee 



ATTTjrj 



(23) 
where 77 is the viscosity of water [56]. We evaluate the 
coe fficien ts of A^^'^B"^ from these expressions. Sec- 



tion 



II D 



provide a similar estimate of A'^ (which is consis- 
tent with experiments as shown in Sec. [VC I, and the co- 
efficients of A'^j^^B'^ are found to be at least four orders 
of magnitude larger than those of A'^. We can therefore 
diagonalize M perturbatively in A'^x- Since A^''y^^B^x 
has one vanishing eigenvalue (see Appendix E I , D3 is 
much smaller than Di and D2. We write it and the as- 
sociated eigenvector to lowest order in A'^x- 



D.= 



A 



det(x) 

Xl,l Xl,2 


-, X3 = 


-X-' 


I 


\ 1 


X2,l X2,2 









(24) 



where the denominator of the second term in D3 is the 
(3,3) cofactor of matrix x- 

This slow mode can be interpreted as follows. The or- 
ders of magnitude calculations given above show that 7"^ 
is large, meaning that the friction of water against the 
tube is very weak. Therefore, according to Eqs. (l5| the 
tube quickly (although in hydrodynamic times oc q^"^) 
relaxes to a state of constant tension and torque Vtr = 0, 
Vr = 0. Anticipating on the results of Sec. |VC| we esti- 
mate that this regime is reached in a few tens of millisec- 
onds for a tube of 10 /im. In experimental situations close 
to that of Ref. [5^, this relaxation is much faster even 
than the injection of GTP in the experimental chamber. 
Therefore, when interested in observable time scales, one 



should consider that the two fast mode s of Eq . ( 21 ) are al- 
ways at equilibrium. Using Eqs. (19a I, (19b I, we deduce 



that p and h have the following spatially homogeneous 
values: 



p = o- 
h ^ -' 



--6 



is 



A/i, 



(25a) 
(25b) 



where a and r are independent of z and are fixed by 
the boundary conditions imposed on the tube. Note that 
introducing GTP in the system, thus changing the value 
of A/i, is equivalent to applying a force — ^^A/x and a 
torque — i^eA/i to the tube. 

If the two fast modes are considered at equilibrium, the 
tube dynamics can be described by the evolution equa- 
tion of the projection of the state of the system onto 
the third mode. In our approximation, this projection 



is 6^1 [x- 



'3,3' 



The equations of motion of the system 



therefore reduce to a single diffusion equation whose dif- 
fusion coefficient is the smallest eigenvalue of M: 



dt5<i> = DsV^d^. 



(26) 



IV. SUSCEPTIBILITY MATRICES 
DESCRIBING EXPERIMENTAL SITUATIONS 

Although this is already true in the general case of 
Eq. (211, it appears even more clearly in the simplified 
Eq. (24) that a full understanding of the dynamics re- 
quires an expression of the susceptibility matrix x- Be- 
fore proposing such expressions, we would like to com- 
ment on the nature of the assumptions they imply. Un- 
like in the previous sections, where only well-controlled 
approximations based on orders of magnitude and the 
symmetries of the system are used, the calculation of x 
requires an explicit expression for the free energy of the 
tube. As emphasized in the introduction, such a mi- 
croscopic description is difficult given our limited knowl- 
edge of the mechanics of dynamin. Nevertheless, since 
the models that we develop in this section are equilib- 
rium models of the tube, all the information about the 
non-equilibrium behavior of the system is still being cap- 
tured by the phenomenological coefficients introduced in 
Section |IIC[ This means that we do not make any as- 
sumptions on the microscopic details of the dissipation 
mechanisms. 

In the following, we first define a microscopic 
parametrization of the dynamin/membrane tube. Then 
we propose three equilibrium models of the tube, aimed 
at describing the experimental situations of Refs. [5^125] 
and [131 mi HSl [TBI [23], where different types of lipids 
were used as templates for dynamin assembly. 



A. Microscopic parametrization 

For the sake of simplicity, let us start by idealizing 
the geometry. We first assume that the membrane is 
infinitely thin. It is confined to a roughly cylindrical 
shape by the dynamin helix but small deviations from 
this shape are allowed in the following. The energy cost 
of such deformations is fixed by the membrane's stretch- 
ing and bending moduli ks and kb- The detailed calcula- 
tion of the membrane's bending energy in the geometry 
considered here is presented in Appendix |B] We further- 
more consider the helix as an infinitely thin inextensible 
elastic rod with spontaneous curvature and torsion, such 
that its equilibrium shape is a helix of radius r = vq and 
pitch 2Trp — 2TrarQ. Its elasticity is described as that of a 
classical spring and is parametrized by its curvature and 
torsional rigidities fc„ and kr [SH] ■ All relevant details are 
presented in Appendix [C] 

The assumption of inextensibility of the rod forming 
the helix is the most speculative point of this section. 
Electron micrographs of dynamin helices treated with the 
non-hydrolyzable GTP analogue GMP-PCP suggest that 
the number of dynamin subunits per unit of helix length 
could change upon GTP binding [T5]. We still use the 
inextensibility assumption for simplicity and by lack of a 
satisfactory alternative hypothesis. 

In the remainder of this article and unless otherwise 




Figure 2: Sinusoidal deformations of the membrane out of its 
cylindrical shape are allowed in our model. The wavy black 
line materializes the intersection of the membrane with the 
a; > 0, y = half-plane. The equation of this line is given 
by Eq. ( |27[ ). Note that it is constrained to touch the helix at 
each period of the tube. 



stated, we express all quantities in units of the helix's 
spontaneous radius tq, the mass per unit length po ^^nd 
the typical force needed to stretch the helix /C ~ 2.2 x 
10~*N (see Appendix O. In these units, we define the 
deviations of the radius and pitch of the helix from their 
spontaneous values hy r — I + Sr and p = a(l + Sp). 

Let A be the area per polar head of lipids and Aq its 
value at equilibrium. We define the relative deviation 
of A as a = a ■ Eventually, and although the con- 
finement by the protein imposes an overall cylindrical 
shape on the membrane, we allow it to bend as long as it 
retains its helical symmetry. We parametrize this defor- 
mation by a number u such that the intersection of the 
membrane with the x > 0, y ~ half-plane is the curve 
(see Fig. [2]): 



X = r(l -t- u{cos{pz) — 1)). 



(27) 



Details are given in Appendix [B| 

The elastic properties of the membrane and the he- 
lix are such that when assembled together, they tend to 
deform each other: the equilibrium configuration of the 
tube is different from the spontaneous shapes of the he- 
lix and membrane taken separately. However, we show in 
the following that these are small effects in the sense that 
at equilibrium, (5r ~ (5p ~ a ~ u ~ to a good approx- 
imation. Therefore, at equilibrium, all the mass of the 
tube is concentrated at a radius vq, hence / = por^ = 1. 



B. Rigid membrane model 

In the experiments of Refs. [531 US], dynamin is as- 
sembled on a mixture of non-hydroxylated fatty-acid 



galactoceramides (NFA-GalCer) , phosphatidylcholine, 
cholesterol, and phosphatidyhnositol-4,5-bisphosphate 
(PtdIns(4,5)P2). The proportions of these Upids are such 
that even in the absence of dynamin, they spontaneously 
form nanotubes with a diameter comparable to that ob- 
served for dynamin-coated tubes. Here we investigate 
a suggestion made in Ref. [55], namely that these lipid 
nanotubes are very stiff. Consequently, we model them 
as rigid cylinders {u = 0) of fixed radius and area per po- 
lar head. The last two conditions are imposed by writing 
the free energy of the tube as 



fo = fh + 



Kr^ 



' X 2. 

-or 



^2o 2 

2 ^ ' 



(28) 



where ft is the elastic energy of the helix, as calculated in 
Appendix [C] The assumptions 5r ^ Q.a ~ Q are enforced 
by taking the limit koo,k'^ —>■ oo. Using the expression 
Eq. ( |C4| for fh and in the limit fcoo, k'^ -^ oo, one induces 
no change in the dynamics by replacing Eq. ( 28 1 by 



foiSr, Sp, a) ^ 'fSr' + 'fSp' + 'fa'. (29) 

Since the rod forming the helix is inextensible, its mass 
density is proportional to the rod length per unit length 
of the tube. The membrane mass density, on the other 
hand, is proportional to the radius of the membrane 
cylinder and inversely proportional to the stretching rate 
of the membrane: 



Ph = $0 



^\/r^ + 



pz 



pV^ 



(1 - $o) 



r 



1 + a 



(30a) 
(30b) 



Combining these equations, one obtains expressions for 
Sp and S^. One also notices that Uze = 1/p— 1/a, hence 
the first-order expressions: 



Sp 



1 — ^ Sr 



1 + a^ 



1 + a^ 



v-zd = 



Sp 



(5$ = (f>o(l - <f>o) 



o?Sr 



2 (1 - $o)a (31a) 

(31b) 

2 .r^+«V(31e) 



1 + a^ 1 + a^ 



Combining Eqs. (29) and (31 1, one obtains the func- 
tion fo{Sp,UzeiS^). The susceptibility matrix is essen- 
tially the matrix of second derivatives of this function 
(see Eqs. (A3 1 for a more rigorous statement). Taking 
the limit k^ , fc' 



i- CX3, we finally find: 



lim 



Arm 



^pp 



(l+a2)2 

fg(l-fo) 
(l+a2)2 






fo(l-fo) 

(1-|-q2)2 

■t'o(l-'^o) 
a(l+Q2) 

M(l-*o)' 



' u V -■- ^ u ; ^ (J V ^ ^ yi/ 



(32) 



Soft membrane models 



and the volume of water enclosed by the tube: 



In many in vitro experiments, dynamin is assembled on 
lipid bilayers containing no cholesterol and which spon- 
taneously form lamellar phases or vesicles in the absence 
of dynamin [T31 HH [T5J ^E\ [53] . From these two observa- 
tions, we can presume that they are much softer than the 
lipids studied above and that their spontaneous curvature 
is zero or at least negligible compared to the curvature 
imposed by the dynamin coat (~ 10* m~^). For these 
lipids, the microscopic variables 5r, Sp, a and u can there- 
fore all take non-zero values. The free energy of the tube 
is therefore the sum of three terms: the spring elastic en- 



ergy of Eq. ( C4 1 , a simple quadratic membrane stretch- 



ing energy with stretching constant ke and the membrane 
bending energy of Eq. (B8l. To second order: 



5r^ 
/o — {Kr + 27rfcb) — - -I- krpSr Sp + k. 



Sp"" 



,2 



+2Trkg h TTki,(—Sru) + fc„ 

+Trkh(—Sr + u). 



(33) 



Minimizing this free energy with respect to Sr, Sp and u, 
we find that at equilibrium 



Spcq 



kb 
JCtq 



-^cq 



a 



2 X 10"'' < 1 



2 X IQ-^ < 1, 



(34) 



(35) 



h h 

f^rpj '^pp 



where we have made use of the fact that kr 
kf^, kr ^ /C in dimensionless units (see Eqs. ( |C3[ )). We 
have considered a typical bending modulus kb ~ 10 ksT 
[39J and estimated a ~ 0.2 from Ref. [23". These orders 
of magnitude show that the linear terms in /o are very 
small and therefore we neglect the last term of Eq. ( 33 ) in 
the following. In other words, we use the approximation 
that at equilibrium the spring assumes its spontaneous 
radius and pitch Toq = 1 , Poq = ol and that the membrane 
is an unstretched cylinder of radius 1 (since Ocq = from 
Eq. (|33l) and Weq = 0). 



As above, we want to express /o (now a function of the 
microscopic variables Sr, Sp, a and u) as a function of the 
hydrodynamic variables Sp, u^g, S^. Since there are four 
microscopic and three hydrodynamic variables, finding 
a unique relationship between the two sets of variables 
seems impossible at first sight. However, there exist con- 
straints on the microscopic variables that we have not 
yet been expressed. To understand these constraints, let 
us calculate two quantities to first order with the help of 
Eqs. (B5l and (B6): the mass density of lipids, which is 



the ratio of the surface area covered by the lipids to their 
area per unit mass: 



Plipids 



= 27r(l + Sr — a — u) 



(36) 



V^TT{l + 2Sr-2u). 



(37) 



If all four microscopic variables were independent, pupids 
and D would be independent as well. However, since the 
membrane tube is filled with water, allowing for a change 
of at constant piipids implies a fiow of the water inside 
the tube relative to the lipids. We estimate the typical 
time scale associated to this fiow to be that of a Poiseuille 
fiow inside a tube of radius tq driven by a pressure dif- 
ference /C/rg and over a distance L = 10 /um: 



tp, 



oiscuillc 



877(7rrgL) L 



m^is. 



(38) 



Therefore, on time scales t <ti ipoiscuiiic, the relative fiow 
of membrane and inner water is insignificant. Conse- 
quently, the ratio of mass density of membrane to mass 
density of inner water /9iipids/(0pH2o) has to be a con- 
stant. Conversely, on time scales t ^ ipoiseuiiie, one can 
consider that the fiow of water inside the tube has re- 
laxed, hence /Oiipids and D are independent variables. On 
time scales t ~ ipoisouiiio, the situation is more complex 
and a correct hydrodynamic theory would involve not 
two, but three different fiuids: the helix, the membrane 
and the inner water. Such a treatment would obviously 
be quite heavy and relevant only on experimentally unob- 
servable time scales. In the following, we therefore only 
calculate x in the two limiting cases t <C ipoiseuiiie and 

t ^ t^Poiscuillc- 



1. Short time scales: t ^ tpoiseuuie 

In this limit, no relative fiow of membrane and inner 
water is possible and the ratio piipids/(f/0H2o) is a con- 
stant (PH2O — lO'^kg.m^"^ - the mass per unit volume of 
water - is considered a constant). Using Eqs. (36 1 and 
this yields to first order: 



(371 



Sr + a — u = 0. 



(39) 



On top of this constraint, one can write three equations 
relating the microscopic variables to the hydrodynamic 
variables. Since the inner water and membrane cannot 
fiow relative to each other, we treat them as a single fiuid, 
which we label "fiuid to", hence /?,„ — pupids + OpHap- 



Similarly to Eqs. (30 1 and using Eqs. (36 1, (37) and (39 1 



one can write to first order 



Ph = *o 



a\/r^ -\- p^ 
p^/1 + a^ 



{\ - ^q){\ + Sr ~ a - u). 



(40a) 
(40b) 
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Moreover, one still has u^e = 1/p— l/a, hence up to first 
order 



5p = 1 - 



5r — 



1 + a^ J ' 1 + a^ 

(1 -$o)a-(l -*o)w 

dp 

a 



S<P = $o(l-$o)(- 



dp 



1 



1 



(41a) 
(41b) 

+ a + u] .(41c) 



Combining these and Eq. ( 39 1 yields a unique relation be 



tween the microscopic and hydrodynamic variables. We 



differentiate the free energy of Eq. ( 33 1 as a function of 
the latter, yielding an expression for Xsm- More details 
are given in Appendix [P] 



2. Long time scales: t ^ tpoiseuuie 

In this limit, the flow of water inside the tube has al- 
ready relaxed and therefore /Oupids and are independent 
variables. Consistently with the hydrodynamic approach 
used in this paper, we consider that the microscopic state 
of the system has the lowest free energy compatible with 
the values of the hydrodynamic variables. This yields the 
following constraint: 



dl 

du 



Q -^ u 



5p, u^Q , (5<& 



Trfcfe 



6r - 



27rfc, 



-a. 



(42) 



As above, this constraint yields a unique relation between 
the hydrodynamic and microscopic variables. Fluid m 
now repre sents only the membrane: p,„ = Pupids- How- 



ever Eqs. (40b I and (41) remain valid. As above, xi 



t» 



obtained by combining them with the constraint Eq. ( 42 ) 
and the second derivatives of /o . See Appendix |D for 
more details. 



V. COMPARISON TO EXPERIMENTAL 
RESULTS 

In this section we use electron microscopy data to eval- 
uate the active force and torque generated by the tube 
when supplied with GTP. Using these results, we show 
that the change of conformation of dynamin is expected 
to depend strongly on whether it is assembled on a soft 
or rigid membrane tube, which could account for seem- 
ingly contradictory experimental results. We then turn 
to the tube dynamics and show that although currently 
available experimental data do not allow a detailed com- 
parison with our theory, the time scales involved are in 
agreement with our predictions. 

The numerical estimates of this section are based on 
the typical values h ~ 10 ksT ~ 4 x 10"^° J [3^ and 
kg ~ 0.25 N.m~^ |35j. ry^atcr = 9 x 10~''Pa.s, and mea- 
surements show To ~ 10 nm, a ~ 0.2, r^ — 25 nm f23| and 



/C ~ 2.2 X 10"^ N (see Sec.|C3|. The water friction and 
helix elastic constants are calculated from Eqs. ( pS] ) and 
(ICSf. We also use po ^ 7.5 x lO^^kg.m^i and $o ^ 0.5 



(see Sec. II D). 



A. Active terms 

According to the symmetry arguments developed in 
Sec. |IIE| exposing the tube to GTP yields the same 
deformation as applying a force — ^^A/i and a torque 
— ^gA/i to it. Making an analogy with a spring submitted 
to a force and torque, we expect a uniform change of ra- 
dius and pitch for a dynamin helix incubated with GTP 
for a very long time. Ignoring fluctuations, this is consis- 
tent with experimental data [T3 | [M l [T5l ITB t [23 l [24l [25] . 

Let us first turn to Ref. [23], where the GTP analogue 
GMP-PCP is used. As discussed is Sec. |IVC| one should 
describe this system with x^^ on long time scales. In 
those experiments, the changes of pitch and radius of the 
dynamin helix are measured to be 



Ar — lim Sr : 

t — >oo 



-0.5, 



Ap — lim Sp ' 

t — >oo 



-0.31, (43) 



which is compatible with the results of Refs. [131 [UJ [T5J 
[TB] . Knowing Xsm> ^'^ ^^^ Ap, we deduce the ampli- 
tude of the active terms ^zAjj, and £,eAfi, just like the 
force and torque exerted on a spring can be deduced 
from its elastic moduli and the amplitude of its defor- 
mation. Considering that no external force or torque are 
exerted on the tube (cr = 0, r = 0) and assuming that the 
dynamin-covered portions of the tube are in contact with 
other regions with which they can freely exchange helix 
or membrane such that fie = 0, we combine Eqs. (20 1 
and ( 25 ) to find 




V Asm ) 







(44) 



Moreover, according to Appendix |D] the left-hand side 
of this equation is a known function of Ar, Ap and Aa — 
limt^oo a. We solve Eq. (44 1 in £,zAfi, ■feA/i and Aa and 
obtain 



C^A^~ 5.8x10"^ N, ^eAfi' 



1.3x10"^* N.m, Aa~0.39. 
(45) 



B. Variability of dynamin's conformation after 
GTP hydrolysis 

In contrast with the results presented above, Refs. jUJ 
?I5\ report that the radius of the dynamin helix does not 
change upon incubation with GTP and that its pitch does 
not decrease but increases, yielding Ar' ~ 0, Ap' ~ 0.7 
(in the following, the dashes denote the deformations as- 
sociated with these references). It is possible to account 
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for those apparently contradictory results without resort- 
ing to biochemical arguments (differences in the type of 
nucleotide or dynamin used...) but only by the mechan- 
ical properties of the lipids. 

We assume that the equilibrium properties of the tubes 
used in these experiments are well described by Xrm, as 



f«T 



T. « t « T, X, « t « T, 



discussed in Sec. IV B Assuming that the biochemistry 
of the tubes considered here is the same as in Refs. [13, 
[HI IT51 \W[ [23] implies that the active terms have the 
values given in E g. (|45| ). Again assuming that /ig = 0, 
we combine Eqs. 



20rand (25 1 to find 




(Xrm)" 







(46) 



It is clear from the assumptions of Sec. |IVB[ that com- 
bining this equation with Eqs. (31 1 and (32) yields Ar' ~ 
Aa' — 0. More interestingly, 



Ap' = - 



auzo 






Ap 



2{1 - a'){k, - kr) 
Ak^a^ + kr{l - a^)^ 



Ar. 



(47) 



Therefore, we predict that the pitch increases {Ap' > 0) 
if and only if 



/Cf,- 



> 



2(1 ~a'^)Ar + Aa'^Ap 
2(l-a2)Ar-(l-a2)2Ap 



1.5. 



(48) 



This condition is not satisfied by the cylindrical rod 
model leading to Eq. ( C3 ) . One should however temper 



this result by considering the crudeness of this model, the 
limited applicability of our small deformation formalism 
to the high nucleotide concentration experiments consid- 
ered in this section, as well as the rather large uncertainty 
on several numerical values used here. We therefore con- 



sider Eq. ( 47 1 as a proof of principle that a shrinkage of 



radius on a soft membrane is compatible with an increase 
of radius on a rigid membrane. 



Time scales 



Turning to the dynamics of the tube as described in 



Ref. [26 , we apply the perturbation scheme of Sec. Ill 
to Xsm using the typical numerical values presented 
throughout this paper. We furthermore assume that the 
size and boundary conditions of the system are such that 
the smallest wave vector allowed is gmin = 27r/(60/im) 
|40j . Eq. (21 ) implies that the deformations characterized 
by Q'min dominate the long-time relaxation of each of the 
three hydrodynamic modes of the system, yielding three 
relaxation times r^ — 2t: / Diq^^^, i = 1,2,3. For simplic- 
ity and without loss of generality, we discuss only these 
deformations in the following. Finally, we assume that 



,t« 



Asm y i --° ** 




^<<t 






t» 



t:,= 120|j,s^ JT2=37ms 'j T3= 2.6 s 




Figure 3: Illustration of the dynamics of the tube generated 
by the susceptibility matrices Xsm and Xsm- Note that to a 
good approximation, the field of deformation of the tube is 
independent of z during the lag phases between the relaxation 
of the chronologically well-separated hydrodynamic modes. 
The amplitude of the changes of conformation are calculated 
from Eq. (451. The transparency of the membrane illustrate 
its stretching and is proportional to 1 + a. The thick black 
arrows represent the expected changes of conformation based 
on the comparison of the r^s with fpoiscuiiic (see text). 



one end of the tube is in contact with a reservoir imposing 
the boundary condition p,e{z = 0) = 0. The relaxation 
time scales are found to be well-separated: ti ~ 120 /iS, 
T2 ~ 37 ms and t^ ~ 2.6 s. This retrospectively validates 
the perturbative approach of Sec. in) 

Comparing ti to tpoiscuiUo (Eq. ([38])), we find that x*^ 
is probably not a good description of the tube on this 
time scale and that some intermediate matrix between 
Xsm and Xsm should be used. Looking at Fig. [3] however, 
we realize that the dynamics generated by the two matri- 
ces are not very different and that using one or the other 
does not make much difference at our level of description. 
On the other hand, it is clear that the transformations 
characterized by T2 and T3 must be described using Xsm • 

In agreement with Sec. |III| ri and T2 are smaller than 
the time needed to inject GTP in the experimental cham- 
ber (typically a few tenths of seconds) and are therefore 
experimentally unobservable: the approximate Eq. (26 1 
is sufficient to describe current experiments. 



On long time scales, our theory predicts an exponen- 
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tial relaxation of the helix with a longest relaxation time 
Ta ex l/(xsmAg^jin) (see Eq. (24 1). This behavior was 
indeed reported in Ref. [25] anothe value T3 — 2.6 s pre- 
dicted by us matches the measured value within exper- 
imental uncertainty. Note however that our choice of 
<Zmin is to a Certain extent arbitrary. Therefore, further 
measurements are required to obtain a precise value of 
A. Still, even only an order-of-magnitude agreement be- 
tween the experimental and a priori determined value of 
T3 suggests that both the picture of the long-time dy- 
namics of the tube as dominated by the friction between 
helix and membrane (Sec. II D I and our description of 
Xsm are essentially correct. 



VI. DISCUSSION 

In this article we describe the dynamics of long 
dynamin-coated membrane nanotubes typically used in 
in vitro, cell-free experiments. This work is therefore rel- 
evant to the biological membrane severing function of dy- 
namin insofar as we assume that the tube breaking mech- 
anisms are similar in those two cases. This last section 
summarizes and discusses our results in this perspective. 

Our formalism describes several previously unac- 
counted for experimental results. Concerning the statics 
of dynamin, we suggest an explanation for the variability 
in change of conformation obtained by different experi- 
mental groups. Moreover, it is reported in Refs. [251125] 
that long tubes incubated with GTP tend to form plec- 
tonemic supcrcoils, which is consistent with our theory. 
Indeed, in our description, tubes held fixed at both ends 
and provided with GTP are analogous to rods with per- 
sistence length £p under a torque —^eAfi and a com- 
pressive force — ^zA/x. Therefore, they supercoil if their 
length is much longer than the critical buckling length 
-^/fcsT^p/^zA/i ^ a few nm, which is always the case in 
practice. Moving on to the dynamics of the tube, we 
show that the longest-lived and only experimentally ob- 
servable internal relaxation phenomenon of the tube is an 
effective friction mechanism between dynamin helix and 
lipids. From this we conclude that the internal dynam- 
ics of the tube can be approximated by a single diffusion 
equation, which accounts for the exponential relaxation 
observed in Ref. [25] . 

This very simple form for the dynamics of the tube im- 
plies robust features that could be tested experimentally: 
first, the dominant relaxation time should scale like the 
square of the length of the tube; second, at long times, 
the rotation frequency of the tube should have a sinu- 
soidal dependence in z, which is the shape of the slow- 
est eigenmode of the diffusion equation. Investigating 
this latter prediction will provide insight into the bound- 
ary conditions on the tubes, which could prove useful for 
making further predictions and may be relevant to tube 
fission. 

In our formalism, nucleotide addition is equivalent to 
exerting a force and a torque on the tube. Hence, af- 



ter relaxation of the transient regime, incubation with 
GTP should induce a homogeneous strain of the dynamin 
helix. Therefore, two points separated by a distance z 
should rotate relative to each other a number of times 
proportional to z. This was confirmed by preliminary 
experimental data |40j . This result is intimately linked 
to the assumption of non-polarity of the dynamin helix, 
which suggests that the subunits of the dynamin helix 
themselves are apolar, as proposed in Refs. [T^ITB] . 

We now discuss the assumptions used to derive our 
formalism. The most important of these is our use of 
the large-system limit L » tq. As discussed in Sec. [IT] 
it is obviously correct when applied to the in vitro, cell- 
free experiments considered in this work. Unfortunately, 
dynamin collars observed in vivo are much shorter - typi- 
cally two to three helical repeats [B]. However, we believe 
that our concepts of friction between helix and membrane 
and GTP-induccd force and torque can be readily trans- 
posed to short tubes. One could also be concerned that 
on small length scales, non-hydrodynamic relaxation phe- 
nomena occur on the same time scales as the relaxation 
phenomena we discuss here and therefore interfere with 
our picture of the relaxation of the tube. To address this 
point, let us note that according to Ref. [25], breaking 
long (~ /im) tubes takes seconds, which is much longer 
than any reasonable non-hydrodynamic relaxation time 
for this system. Equivalently, we can say that the tube 
does not break in short (non-hydrodynamic) times, from 
which we conclude that non-hydrodynamic internal re- 
laxation phenomena are not essential to tube breaking. 
Therefore, although the in vivo situation is undoubtedly 
more complex than that considered here, we argue that 
our description of the internal dynamics of the tube is 
sufficient to study its tube-severing function. 

By writing constitutive equations for the tube, we as- 
sumed it to be weakly out of equilibrium. Concerning 
the friction and viscosity-related phenomenological coef- 
ficients (those of Eqs. (17l), experience shows that this 



requirement is not very stringent [28 , and our consti- 
tutive equations are likely to give a good description of 
the system in most situations. Chemical systems, on the 
other hand, typically operate far from equilibrium. In 



this regime, writing the active forces of Sec. HE as A/i 
does not yield the correct dependence on the concentra- 
tion of GTP. Better results would probably be obtained 
by using instead 1 — q-^^^/>^bT ^ which is characteristic of 
molecular motors |41j . Forthcoming experiments involv- 
ing low levels of GTP are expected to better fall in the 
domain of applicability of our theory. 

Assumptions of small deviations of the tube from its 
initial state are also used when deriving the suscepti- 
bility matrices x- Again, these are formally correct at 
small concentrations of GTP. However, given the fact 
that those matrices involve uncertain microscopic as- 
sumptions, we regard them as nothing more than reason- 
able examples anyway. Therefore, we do not expect them 
to yield quantitative results. More reliable information 
about the characteristics of these matrices could be ex- 
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tracted from micromechanical experiments on dynamin- 
coated nanotubes. 

We now turn to dynamin-induced tube breaking mod- 
els from the literature. We do not discuss the purely 
biochemical model of dynamin as a regulatory GTPase 
- which has consistently been regarded as unlikely over 
the past few years [TJ [21 [22] - and rather concentrate on 
the mechanochemical models. Depending on authors, the 
critical feature leading to tube breaking by dynamin has 
alternately been proposed to be a change of radius [TO] , 
of pitch [24 , or membrane bending [42 . As can be seen 
from Sec. jIVCj all of these deformations fit very natu- 
rally in our theoretical framework. Moreover, we showed 



in Eq. ( 45 1 that the change of conformation of dynamin 



typically induces a significant stretching rate Aa of the 
membrane. We would therefore like to attract attention 
on this fourth type of deformation of the tube, which 
might play an important role in tube breaking. 

Unlike previous models, this work does not rely on 
detailed assumptions about the tube's conformational 
changes. Instead, we predict them by optimally exploit- 
ing the experimental data in the light of thermodynamic 
considerations, conservation laws and symmetry argu- 
ments. 

Several models have also been proposed for the cou- 
pling of GTP to dynamin activity: GTP hydrolysis could 
induce a concerted conformational change |25 , while 
some results suggest that the crucial step is the binding 
of GTP to dynamin |33] and others point to a ratchet-like 
mechanism for its constriction |44j . Since in our frame- 
work the coupling of the energy source to the dynamics is 
deduced from symmetry considerations only, our theory 
is equally valid in each of these cases. 

In addition to including most features previously dis- 
cussed in the literature, our formalism yields novel quan- 
titative insight into the mechanism of tube breaking. It 
would be interesting to further discuss Refs. [HJ H5] . 
where it is assumed that lipids cannot fiow through 
the dynamin-coated region of the tube, in parallel with 
Ref. |1B], where this flow is on the contrary considered 
instantaneous, in the light of our new knowledge of the 
helix/membrane friction coefficient A. 

Furthermore, our hydrodynamic point of view could 
account for some discrepancies between existing models 
and experiments. Adopting the classification of Ref. [2^ . 
we consider models of the "Garrote" class - where the 
reduction of the dynamin radius pinches the membrane 
tube to its breaking - and of the "Rigid helix/Elastic 
membrane" class - where the tube breaks because its 
walls are brought together by a sudden deformation of 
the helix. If taken at face value and applied to a long 
tube, these local models predict a uniform density of tube 
breaks, since what is expected to happen at the neck of 
a clathrin-coated endocytosis vesicle should happen at 
every point of the long helix. Experimentally, however, 
no breaking is observed in such tubes unless their ends 
are firmly attached to a fixed substrate |22j- Moreover, 
attached tubes are observed to straighten upon GTP in- 



jection and then break not at several but at a single point 
|26j . This sensitivity to distant boundary conditions and 
spatially inhomogeneous behavior of the tube motivate 
our description of long-range interactions mediated by 
tube elasticity and the z-dependence of the tube defor- 
mation, which could account for the existence of a pre- 
ferred point of breaking. 

In the "Spring" model [24] as well as in Ref. |45j . 
breaking only occurs at the interface between a dynamin- 
coated and a bare region of the membrane nanotube. 
We can imagine that in long tubes, such defects in the 
dynamin coat either appear during polymerization or 
that the initially homogeneous dynamin coat breaks upon 
GTP injection. It is however very unlikely [30] that the 
dynamin coats of Ref. [2H^ have systematically exactly 
one defect, which would account for the fact that they 
break at most once. Instead, the tube probably often 
starts with either many or no defects. In the former 
case, we have to account for the fact that only one of the 
defects evolves into a full breaking of the tube. In the 
latter case, we must explain the creation of a defect in 
the dynamin coat. It is undoubtedly important to con- 
sider the space dependence of the stresses in the tube to 
answer either of these questions. A mechanism of lipid 
phase separation similar to that of Ref. [37] could also 
help a defect evolve into a full tube break. Indeed, dy- 
namin is known to strongly bind PtdIns(4,5)P2 [48j, and 
could therefore deplete the bare membrane regions in this 
lipid. Also, depending on the flow of membrane through 
the dynamin-coated regions of the tube and on the dy- 
namics of their change of conformation, a bare region 
might be under more or less stress and therefore break 
more or less easily. 

In conclusion, we developed a complete theoretical 
framework suited for the analysis of the statics and dy- 
namics of long dynamin-coated membrane nanotubes. 
We make several predictions concerning the space and 
time dependence of forces, torques, membrane tension, 
membrane stretching and helix deformations. We hope 
that our theory will facilitate the interpretation of forth- 
coming experimental results and help generate and quan- 
titatively test novel hypotheses on the biological mode of 
action of dynamin. 
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Appendix A: DERIVATION OF THE 

HYDRODYNAMIC MODES AND 

PERTURBATION THEORY 

This appendix contains details about the derivation 



(Sec. II F I and simplification (Sec. Ill I of the hydrody- 
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namic equations for the tube, 
vation equations Eqs. (|3|, Q, 
equations Eqs. ( 19 ) we write the hnearized equations of 



Starting from the conser- 
( 10 1 and the constitutive 



motion for the system: 



= A''V 




^dy2 
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h 



(Al) 



B'' + j 
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A" V 
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B^V^ 



(A2) 



where the superscripts "r" and "d" denote matrices of re- 
active and dissipative couphngs respectively. These ma- 
trices read 




A^= -i, "r ,A' = 



-I -i 
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Let us consider Eqs. (Ill, (12 1 and (13). They imply 
that 



P 



d{fo/p) 



dp 



h = 



dfo 



*,Mi8 



du. 



Pe = 



p,<E> 






P-M-^e 

(A3) 
where the temperature T is considered constant. The 
fields p, h, fie are therefore functions of Sp, Uze, and (5$ 
only. Close to equilibrium this dependence can be lin- 
earized, yielding the definition Eq. ( 20 1 of the suscepti- 
bility matrix x- Using this definition, we can now obtain 
a closed equation for the hydrodynamic variables. In 
the presence of friction against water (i.e. if 7 is posi- 



tive definite) , the left-hand side of Eq. ( A2 1 is irrelevant 
in the hydrodynamic limit, as shown when going from 
Eqs. (Ill) to Eqs. (Isl. To leading order in the wave vec- 
tor q, this, together with Eqs. (All and (20) yields the 



hydrodynamic modes equation Eq. (21), with 



A'' 




(A4) 



We now turn to the perturbative diagonalization of 
matrix M — (A''j^^B^' + A'^) x, a slight generalization 
of first-order quantum mechanical perturbation theory 
to non-hermitian matrices f49J. The important point is 



that the unperturbed matrix A'^'y~^B'^x has one vanish- 
ing eigenvalue D^ — 0. Indeed, the definitions of A^ and 
B^ imply that 



^'^7-15'- = 




(A5) 



where the question marks stand for non-zero coefficients. 
Clearly, the vector x^ of Eq. ( 24 ) is an eigenvector of M 
associated with D^. To lowest order in A'^x, D^ is given 
by Eq. (24 1 and x^ is the associated eigenvector. 



Appendix B: MEMBRANE GEOMETRY AND 
BENDING ENERGY 



The calculations of this section are inspired by those of 
Ref . [4^ . In this appendix, we calculate the bending en- 
ergy of an infinitely thin membrane with no spontaneous 
curvature surrounded by a helical scaffold of radius r and 
pitch 2TTp. We assume that this scaffold imposes two con- 
straints on the membrane: first, the membrane has the 
same helical symmetry as the scaffold; and second, the 
membrane is attached to the scaffold and must there- 
fore touch it at every point. Under these constraints, the 
membrane radius as a function of the angle 9 € [—00, 00] 
of cylindrical coordinates and the elevation C G [0,27rp] 
from the scaffold reads (see Fig. B : 



r™(0,C)-r(l + e(C)), 



(Bl) 



where e(0) = e(27rp) = 0. It can be shown that in the e ^ 
1 regime, approximating e by its first Fourier component 
changes all results presented in this article by less than 
1%. We therefore use this approximation throughout: 

r„(0,C)=r(l + u(cos(K)-l)) (B2) 

The bending free energy of the membrane reads 

'h. 



Frr 



W d6, 



(B3) 



where TL is the total curvature of the membrane and fc^ its 
bending modulus. The integration runs over the surface 
of the membrane. For the configurations considered in 
this paper, the dependence of fcf, on the area per polar 
head of lipids can be neglected [50j. We now calculate 
Fm as a function of r, p and e using differential geometry 
[ST]. To second order in e, the surface element of the 
membrane reads 



dS = 1 + e 



1 



p- 



(9^e)^ rd6ldC. (B4) 



Integrating this surface element and using Eq. ( B2 1 , we 
calculate s, the membrane area per unit of z- length of 
the tube to first order in u: 



27rr (1 — u) . 



(B5) 
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torsion f35^. Its elastic energy therefore reads 



Figure 4: Parametrization of a surface confined by a he- 
lical scaffold. Any point whose cylindrical coordinates 
{rm,Om,Zm) cau bc written as Vm = r{l + e(C)), dm = 6, 
z-m = p6 + C, with 9 G [—00,00], ( G [0,27rp] belongs to the 
membrane. 



Similarly, to first order in u the volume enclosed by the 
membrane by unit of z-length is 



Fh = 



(^'^inii)-Kor+'^irii)-ror^d£ 



yjr"^ + p^ I ^K f r 



P 



Tq 



fcr / P arp 

2 \r2+p2 Tq + (aro)2 



1 1 d2, (CI) 



= Trr^ (1 - 2u) 



(B6) 



where the first integral is calculated over the curvilin- 
ear length of the rod and the second one over the z- 
coordinate (see Fig. Ilia)). In the second expression we 
replace the curvature k(^) and torsion t{(,) by their values 
for a spring of radius r and pitch 2'Kp. We also choose 
the spontaneous curvature kq and torsion tq such that 
the ground state of the rod is a helical spring of radius 
tq and pitch 2i:arQ. 

The dynamin helix binds to the membrane through a 
specific domain, the PH domain. Let u{t) be the unitary 
vector field defined on the helix that always points in 
the direction of the PH domain. Since the PH domain 
always faces the membrane, u(£) always faces the inside 
of the helix. Hence \i{t) = N(£), the normal to the helix. 
Consequently, according to Ref. [S^, the twist density of 
the helix is exactly equal to its torsion, which allows us 
to write Fh as a function of n{t) and t(^) only. 



To second order in e. The total curvature of the mem- 
brane reads 



n = 



-'(l-^)+fl + 5)K9^) 



9 9 



r r 2 V r"^ 



2^-e{die) + 1 - ^ r{d^eY. (B7) 



Performing the integration of Eq. (B3) using Eqs. (B2| 



( B4 1 and ( B7 1 , we find the bending energy per unit of 



z-length of the tube 



Jm 



3 3a-2 



u 



^(i + «+li + ^ + "-',, 



- 7rfc6(l + u) + kuu^ 
r \ 2 



(B8) 



This last equality defines k^ 



1. Curvature and torsion coefficients of a rod. 

In order to calculate the curvature and torsion mod- 
uli of the rod, we consider it as a rod of cross-section 



TT I 



-rp - 



where re is the outer radius of the dynamin 
coat. We first define the typical force needed to deform 
the helix 



/C = 



TTEi 



zTsl 



(C2) 



where E is the Young modulus of the rod. Assuming 
that the Poisson ratio of the rod is 1/2, its curvature and 
elastic moduli read 1531 
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(C3a) 



(C3b) 



where the equalities on the right are valid if all quantities 
are expressed in units of ro, pn, IC. We use these scaled 
units in the rest of this section. 



Appendix C: ELASTIC PROPERTIES OF THE 
HELIX 



2. Elastic energy of a spring 



We describe the elasticity of the dynamin helix as that 
of a simple rod with constant spontaneous curvature and 



Writing r — I + 6r and p = a(l -I- dp), the integrand of 



Eq. (CI I (i.e. the elastic energy per unit z-length) reads 
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to second order in Sr, Sp 

Jh 



(l + a2)7/2 

+2 (1 - a^) {k^ - k-r)SrSp 
+ (4a2fc^ + (1 _ a^kr 



Sp-" 



2 



6p^ 



f^rr "^ ~r f^fpOT OJ) ~t f^pp jr~ 



(C4) 



This equality defines krr, krp and kp 



3. Persistence length of a helix and experimental 
determination of IC 

We now consider the possibihty for the central axis of 
the helix to bend with a radius of curvature TZ^ tq. For 
simplicity, we assume that the radius of the helix remains 
constant under this deformation. The pitch is allowed to 
vary over one period of the helix to minimize the energy 
of the system. The persistence length of the helix can be 
defined by expanding its elastic energy in powers of TZ: 



Fh{n)-Fh{n = Q)^kBT j ^dz 



ngEq. (CI), this yields: 
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Appendix D: SOFT MEMBRANE 
SUSCEPTIBILITY MATRICES 



For clarity's sake, we regroup some expressions associ- 
ated with the models of Sec. II V CI here. In both limits 



discussed in Sec. IV C (short and long time), one com- 
bines the relations Eqs. (41 1 between the three hydrody- 



namic and four microscopic variables with a constraint 
(Eqs. (39 1 and (42) respectively) to find a unique lin- 
ear relation between the hydrodynamic and microscopic 
variables, which we write: 




Qs 




(Dl) 



In a unit system such that po = 1, Eqs. (|20[) and (A3 1 



imply that x is the matrix of second derivatives of the free 
energy of Eq. (33) as a function of {Sp,Uze,S^)- Hence 



X«m= ((Qsm) ')^ifsm(gsm)"'. (D2) 

The expressions for the matrices implicated in this for- 
mula depend on the time scale considered. They read 



f^rr ' ^uu ^rp '^uu '' "-b 
jV^„, ^ I >^rp '^pp U 

kuu T^kh U Zirkg + ky^u 



(D3) 



In Ref. [54], the persistence length of a nan- 
otube of brain polar lipids and phosphatidylinositol-4,5- 
bisphosphate coated with rat brain dynamin is measured 
to be Ip = 37±4/xm. From this and Eq. (C6l, we calcu- 
late 
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Using Eq. ( C2 1 , we estimate E ~ 220 MPa, a value of the 



same order of magnitude as the typical Young modulus 
of proteins E = 2GPa [SS]. 
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Neither xim ^'^^ Xsm have compact explicit expressions. 
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Note that as far as Eq. (261 is concerned, the exact values 
of the coefficients of 7^ do not matter as long as they 
are large. Therefore, there would be no point in trying to 
improve Eq. (23 1. 



